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Abstract 

The interplay between disorder and interactions is a leit-motiv of condensed matter 
physics, since it constitutes the driving mechanism of the metal-insulator transition. 
Bose-Einstein condensates in optical potentials are proving to be powerful tools to 
quantum simulate disordered systems. We will review the main experimental and 
theoretical results achieved in the last few years in this rapidly developing field. 
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1 Introduction 



In Nature many processes occur in an ordered way. Indeed, ordered config- 
urations are often the ones minimizing the total energy of the system. A 
prominent example of this tendency towards order is given by the growth of 
a crystal, where the atoms arrange themselves in a spatially periodic configu- 
ration building up an ordered lattice. The physics of transport of electrons in 
a metal heavily relies on the periodicity of this lattice. However, when crys- 
talline solids are studied on a sufficiently small length scale, one realizes that 
impurities and defects are always present, which may affect in a substantial 
way the transport of the electrons. Disorder is indeed an intrinsic property of 
all the real systems. In the last 50 years the effects of disorder on transport 
phenomena have been extensively studied in the context of both statistical 
and condensed matter physics. 

Despite the very general interest in understanding the physics of disorder in 
condensed matter systems, still many questions remain open and unsolved, 
even from the theoretical point of view. As a matter of fact, the theoretical 
description of periodic systems, as perfect crystals, is much easier than the 
one for disordered system as disordered lattices or glasses. The problems that 
arise are related to the fact that the effects of disorder cannot be theoretically 
treated in a perturbative way: even a small amount of disorder can produce 
dramatic changes in the physical properties of the system under investigation. 
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In 1958, P. W. Anderson published a seminal paper (lAndersonl . Il958l ) in which 
he showed under which conditions non-interacting electrons in a disordered 
metal can either move through the system, or be localized. It was soon realized 
that Anderson localization is a much more general phenomenon holding for 
the propagation of generic classical waves in disordered media. Localization is 
a coherent effect that arises from multiple scattering of a wave from randomly- 
distributed impurities and from the resulting destructive interference in the 
direction of propagation. 



Also interactions are well known to induce localization effects, as pointed out 
by N. F. Mott who was able to explain the anomalous insulator behavior of 
some materials when electron-electron interactions were included in the band 
theory. In 1977 P. W. Anderson and N. F. Mott were awarded with the Nobel 
Prize in Physics for their fundamental theoretical in vestigation s of the elec 



tronic structure of magnetic and disordered systems (lAndersonl . 1 19781 : iMott 



19781 ). Following these pioneering works, a strong theoretical effort has been 
devoted in the last decades to investigate the combined role of disorder and 
interactions in the superfluid-insulator transition observe d in many cond ensed- 
matter systems, such as ^He adsorbed on porous media (jCrowell et al. 
thin supercon ducting films (IGoldman and Markovid Il998l ). arrays of Joseph- 



tors (jJiang at al.l . 1 1994 : iBudhani et al 



son j unctions (I Van der Zant et al.l.ll992) and high-temperature superconduc 



1994). 



Ultracold atoms in optical lattices (IMorsch and Oberthaleii . l2006l ) represent an 
extremely powerful tool for engineering simple quantum systems wi th a broad 
tunability of the parameters, thus serving as "quantum simulators" (IFeynmaru . 
19821 ) to reproduce the physics of different systems. The striking advantage 
offered by such atomic systems resides in the unprecedented possibility to 
work with perfectly isolated samples at quasi-zero temperature and to have 
experimental control on most of the Hamiltonian parameters, e.g. the lattice 
depth or the strength of the atom-atom interactions, that can be precisely 
tuned even in real-time. One spectacular demonstration of this opportunity 
has been given by the observation of the superfluid (SF) to Mott insulator 
(MI) transition in a 3D optical lattice (IGreiner et al.l . |2002| ). which pioneered 
the investiga t ion of strongly quantum correlated regimes with ultracold atoms 
flBloch et al.l . l2007h . 



A natural extension of these experiments is the realization of disordered sys- 
tems using ultracold atoms in optical potentials. In this paper we will review 
the recent progresses in this field, that was experimentally initiated in 2004 
with the first investigation of atomic Bose-Einstein condensates in disordered 
potentials. Different possibilities can be followed to produce disordered ultra- 
cold atomic systems. Disordered or q uasi-disordered p otentials can be created 
optically by using speckle patt erns (ILve et al.l. 120051) or multi-chromatic in- 
commensurate optical lattices ( Fallani et al. . 20071 ). These methods allow the 
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production of disordered potentials in which both the spectral properties and 
the amount of disorder are known with very good accuracy and can be easily 
controlled. In addition to the opt ical way, disordered system s could also be cre- 
ated by using a tomic mixtures (Gavish and Castin . 20051) or in homogeneous 



magnetic fields ( Gimperlein et al. . 20051 : Courteille et al. . 2006 ). We will re- 



view these different possibilities together with the illustration of the diverse 
interaction regimes that can be investigated. The first experimental results 
obtained with ultracold bosons in disordered potentials will be presented, dis- 
cussing the state of the art of this newborn field and the perspectives for future 
breakthroughs. 



2 How to produce a disordered potential 



In this section we will present different experimental approaches to the pro- 
duction of disordered potentials for neutral atoms. We will mostly focus on 
two methods allowing the production of complex optical potentials: speckles 
patterns and multichromatic lattices. 



2. 1 Speckle patterns 



The first realization of disordered po t ential s for cold atoms has been obtained 



with speckle patterns (IBoiron et al.l . Il999l ). Speckles are produced whenever 



l ight is reflected by a rough surface or transmitted by a diffusive medium 
( IGoodmanl . l2006l ). We will mostly consider the case of transmission, sketched 
in Fig. la, and we will refer to the scattering device as a diffusive plate. Such 
device can be modeled as made up of many randomly-distributed impurities 
by which the illuminating light is scattered. Since the scattering of laser light 
is mainly a coherent process, the partial waves emerging from the scattering 
interfere and produce a complex distribution of light, called a speckle pat- 
tern, an example of which is shown in Fig. lb. This disordered distribution 
of light can be imaged onto the atoms, producing a disordered potential V{y) 
proportional to the local laser intensity /(r). 



In general, if the wavelength of the light is far detuned from the atomic res- 
onance, no absorption is involved and the resulting mechanical effect can be 
described by a potential energy of the form 



where c is the speed of light, ujq is the frequency of the atomic resonance, F 
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a) 



speckles 



laser beam 



b) 




diffusive plate 



Fig. 1. Production of speckle patterns, a) A laser beam is shone through a diffusive 
plate and the resulting speckle pattern is then imaged onto the BEC. b) Intensity 
distribution of a typical speckle pattern recorded with a CCD camera. 



its radiative linewidth, A = u — uq the detuning, and I (r) the intens i ty dis - 



tribution. This potential is often called dipole potential (iGrimm et al.l . |2000| ) 



It is worth noting that the sign of this potential depends on A, which is the 
only quantity which can take either positive or negative values. In particular, 
when A < {red detuning) V{r) is negative, hence maxima of light intensity 
correspond to potential minima: atoms will move towards higher-intensity re- 
gions. Instead, when A > {blue detuning) V{r) is positive, hence maxima 
of light intensity correspond to potential maxima: atoms will move towards 
lower-intensity regions. 

Speckle patterns represent a valuable way to produce a disordered potential 
in a controlled way. The possibility to accurately measure the statistical and 
correlation properties of the disordered potential comes from the fact that the 
intensity of the speckle pattern can be directly recorded by a CCD camera 
(typically the same one used to image the BEC atoms). In Fig. 2a we show 
the cross section of a typical speckle pattern use d at LENS for the first investi- 
gation of disordered Bose- Einstein condensates (ILye et al.l . 120051 ). Among the 
different quantities characterizing the properties of the speckle field, one can 
define an average speckle height Vs. Different definitions are used in literature, 
however one of the most used in the context of BEC experiments corresponds 
to taking twice the standard deviation of the speckle potentia l V{x) (supposed 
one-dimensional) around its mean value V ( iLye et al.l . l2005l ) : 



Vs 



L/2 



x) — V] dx 



-L/2 



1/2 



An even more important quantity, as we shall see in the following, is the auto- 
correlation length a, giving information on the speckle grain size (iGoodmanl . 
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a) b) 




position [pm] distance [|jm] 



Fig. 2. Production of speckle patterns, a) Cross section of a typical speckle potential 
(the energy difference between the horizontal lines is the average speckle height Vs)- 
b) Autocorrelation integral of the same potential. 



20061 ). This quantity is defined as the rms width of the autocorrelation integral 



G{d) of the speckle potential 



L/2 

G{d)= J V{x)V{x + d)dx^e-^ , (3) 

-L/2 



an example of which is shown in Fig. 2b. The autocorrelation length a de- 
pends on the wavelength of the light, on the nature of the diffusive medium 
producing the speckle pattern and, most importantly, on the optical resolu- 
tion of the lens system used to image the speckle pattern onto the atoms. As 
a matter of fact, the typically determined autocorrelation length is set by the 
diffraction limit spot size of the imaging sy s tem. A detailed description of the 
speckle potential is given in lClement et al.l (120061 ) . where the statistical prop- 
erties of the speckle field are discussed both from a theoretical point of view 
and with the introduction of experimental methods which allow their precise 
determination. 



The random potential produced by a speckle pattern is static. This means 
that the atoms experience just one realization of disorder, which can be repro- 
duced in the same way from one experiment to another. However, shifting the 
position of the diffuser leads to a different realization of the speckle pattern 
that preserves the same spectral and statistical properties. Thus, averages on 
multiple realizations of disorder can be achieved in a simple way. 

Before concluding this section, we note that speckles are intrinsically two- 
dimensional in the plane perpendicular to the propagation axis. Actually, a 
speckle pattern also varies along the direction of propagation of the light. How- 
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ever, the typical correlation size along this direction is much larger. Nonethe- 
less, speckle potentials with different dimensionality can be produced: ID 
speckles can be produced by using cylindrical lenses stretching the speckle 
pattern along one direction, while 3D speckles could be obtained by adding 
speckle patterns coming from different directions. 



2.2 Multichromatic lattices 



As we have seen in the previous section, speckle patterns are a powerful and 
easy-to-implement method to produce random potentials. We have pointed 
out that a crucial parameter of such potentials is the autocorrelation length 
a, which gives an estimate of the minimum length scale below which the po- 
tential loses its random nature and becomes correlated. Typically, this length 
is connected with the diffraction limit dimension at which optical speckles are 
imaged onto the atomic sample. For this reason, the random potential pro- 
duced by speckles is often too coarse-grained (with a of the order of several 
microns), unless one builds a dedicated setup to overcome the usual optical 
access restrictions. Recently, progresses in the realization of speckle potentials 
with autocorrelation len gth below 1 nm have b een achieved in the e xperimen- 
tal groups of A. Aspect fIClement et all . l2006h and iDeMarcol fl2007l ). 



On the other side, having in mind many years of exciting physics with cold 
atoms in optical lattices, we know that optical standing waves can be easily cre- 
ated providing spatial periodicities that can also be smaller than half a micron 
(roughly speaking, one order of magnitude smaller than the autocorrelation 
length of the speckle pattern shown in Fig. 2). This suggests the idea that, 
by combining several optical standing waves with different non-commensurate 
spacings, it is possible to produce complex potentials with very small "grain 
size" . 



The simplest example is given by a bichromatic lattice resulting from the 
addition of two lattices with incommensurate wavelengths. Bose-Einstein con- 
densates i n potentials o f this kind have been first investigated in experiments 
at LENS (IFallani et al.l . 120071 ). where a main lattice with wavelength Ai was 
perturbed by a weaker secondary lattice with wavelength A2, as sketched in 
Fig. 3. The resulting potential can be written in the form 



V{x) = siEri cos^ {kix) + S2ER2 cos^ {k2x) 



(4) 



where ki = 271 /Xi and ^2 = 27r/A2 are the lattice wavenumbers and si and S2 
are adimensional numbers indicating the heights of the two lattices in units of 
the recoil energies Em = h? /{2m\\) and Er2 = ^^/(2mA2), respectively. In 
the limit S2 <^ Si the height of the optical barriers is roughly constant across 
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main lattice V(x) = s^Ej^^COS'ik^x) 

llWWWWWWWWWt 



disordering lattice F{x) = ^jf^^ cos (kjX) 

WWWVWWWWWV^ A-2 = 1076 nrr 

bichromatic lattice F(x) =5|£'j„cos^(A^jX) + 5,£'^2COS^(A:2x) 

= 830 nm 
+ 

X2 = 1076 nm 

Fig. 3. A bichromatic optical lattice. The discrete translational invariance of the 
main lattice is perturbed by the addition of a secondary lattice with incommensurate 
wavelength. 




the whole lattice and it is possible to define a tunneling rate J which only 
depends on the main lattice height si. In this limit, the effect of the secondary 
lattice reduces to an inhomogeneous and non-periodic shift of the potential 
energy at the bottom of the lattice wells (see Sec. 3.2.2 and Sec. 4). 

As we will further discuss in the following, bichromatic incommensurate lat- 
tices are not truly disordered potentials. They differ from both purely random 
potentials and speckle potentials, which exhibit different statistical and cor- 
relation properties. Strictly speaking, they are quasiperiodic potentials, since 
their spectrum is made up of a set of discrete frequencies. However, because 
of the lack of any translational invariance, they can be used to investigate the 
physics of finite-sized disordered systems and study the emergence of quantum 
localization effects, as we shall see in Sec. 3.2.2. 

We note that, since one always deals with finite-sized atomic samples, the 
notion of incommensurability (i.e. the wavelength ratio being an irrational 
number) is a rather sophistic concept, and should be substituted with a more 
practical definition. From an experimental point of view, since the lattice 
wavelengths are known with finite precision, the measurement of the ratio 
A2/A1 always gives a rational number. From a theoretical point of view, it is 
important to consider that the finite size of the systems under investigation 
releases the constraints on the incommensurability: even a periodic potential 
(resulting from a commensurate ratio) does not show any periodicity if the 
system size is smaller than the period. The bichromatic lattice is thus ejfec- 
tively incommensurate provided that the ratio between the wavelengths is far 
from a ratio between simple integer numbers. More precisely, a bichromatic 
lattice can be considered incommensurate whenever the resulting periodicity 
(if any) is larger than the system size. 
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Fig. 4. Disordered potential produced by the interaction of the atoms with a differ- 
ent species. An incoming test particle (white) experiences the collisional potential 
prod uced by randomly distribu ted (black) atoms trapped in an optical lattice. Taken 
from Gavish and Castin ( 20051 ). 



2.3 Other methods 



In solids disorder is often caused by the presence of impurities, i.e. atoms of a 
different kind that randomly occupy the sites of the crystalline lattice where 
atoms of different species were expected. This kind of disorder can be simu- 
lated also in cold gase s by using a mixture of tw o different atomic species, as 
originally proposed in lGavish and CastinI (120051 ). The atoms of one of the two 
species are trapped in the sites of a deep optical lattice. If the filling factor 
(i.e. the average number of atoms per site) is less than unity, only some of 
the sites will be occupied by one atom and the other ones will be empty, as 
schematically shown in Fig. 4. The atoms of the other species, that could be 
weakly affected by the presence of the lattice, feel the collisional interaction 
with the randomly-distributed atoms of the first species. This kind of disor- 
der is spectrally different from both the speckle and bichromatic potentials, 
sin ce it is a binary ki n d of d isorder (yes/no) on top of a periodic backbone. 
In iGavish and CastinI (120051 ) the authors have theoretically investigated the 
possibility to study ID Anderson localization of matt er waves with this sys- 



tem a nd this work has then been extended in 3D in iMassignan and Castin 



( 120061 ). This scheme has not been experimentally realized yet, although first 
experiments with binary ultracold mixtures have been performed (see Sec. 
4.2). 



Anoth er way to introduce disorder in the system has been proposed in lGimperlein et al 
( 120051 ) by using inhomogeneous magnetic fields, e.g. by exploiting the magnetic 
field fluctuations in the proximity of a microtrap caused by imperfections in 



the chip fabrication (IWang et al.l.l2004l) . If the bias magnetic field is kept close 
to a Feshbach resonance (jinouye et al.l . Il998l ). small field fluctuations on top 



of it produce spatial fluctuations in the scattering length characterizing the 
interactions between the atoms. Therefore, this technique allows to introduce 
disorder on the atom-atom interaction strength, rather than on the external 
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potential. In iGimperlein et al.l (120051 ) the phase diagram of interacting bosons 
in the presence of such disorder has been derived, evidencing novel features 
with respect to the phase diagram with disorder in the external potential (that 
will be presented in Sec. 4). 



3 Weakly interacting regime 



One of the most fascinating phenomena characterizing the transport of waves 
in random systems is Anderson localization. This effect takes its name after 
the seminal work of P. W. Anderson in 1958, who identified the fundamental 
role o f disorder i n the metal-insulator transition observed in solid state sys- 
tems ( lAndersonl . Il958l ). Anderson first formulated his localization theory for 
a simple model of particles hopping on a lattice with random on-site energies, 
arguing that above a critical disorder amphtude the quantum states had to 
change from extended to spatially localized. This intuition, together with the 
mathematical tools developed to describe the localization transi tion, led to 
the award of the Nobel Prize in Physics in 1977 (lAndersonl . Il978l ). 



In the following decades, however, it was realized that Anderson localization 
is a much more general phenomenon, holding for propagation of generic linear 
waves in diso rdered media. I n deed, it has been obse r ved f or sound waves and 
light waves (IWiersma et al.l . Il997l : ISchwartz et al.l . 120071 ) . whereas a direct 



observation for matter waves has not yet been possible. 



In the language of wave propagation, Anderson localization arises because 
of interference effects in the scattering of a wave by disordered defects. When 
studying wave propagation in disordered systems, different localization regimes 
can be identified. A precursor effect of Anderson localization is weak localiza- 
tion, which arises from interference effects in multiple scattering events: an 
example of weak localization is given by coherent backscattering, i.e. the en- 
hanced probability of backdiffusion for light incident on a disordered sample, 
owing to the interferen ce between the forward and backward scattering paths 
( IWiersma et al.l . Il995l ) . In the strong scattering limit /c/ ~ 1 (with k wavevec- 
tor and / mean free path between scattering events), these interferences can 
add up to completely halt the waves inside the random medium, resulting in 
strong localization, or Anderson localization. 



Anderson localized states are characterized by the typical exponential decay 
of their tails in the space distribution. In the case of light, this means that the 
intensity is an exponentially decreasing function of the distance travelled in 
the disordered medium. For quantum-mechanical wavefunctions, this means 
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that a localized state \&(a;) can be written as 

^{x)^exp(-^) , (5) 



where ( is the localization length. Generally speaking, the stronger is the dis- 
order the smaller is the localization length. 

In the physics of Anderson localization an important role is played by the 
dimensionality of the system. After the first Anderson's conjecture, scaling 
arguments have been proposed which predict different scenarios with changing 



dimensionality d of the system (lAbrahams et al.l . 119791 ) . For d < 2 all the states 



are localized. For d > 2 a localization transition exists, with a mobility edge 
separating extended states for weak disorder from localized states above a 
critical value. The 2D case is marginal, since the states are localized for any 
amount of disorder as in ID, but the localization length at weak disorder can 
be exponentially large. 

A Bose-Einstein condensate is characterized by long-range coherence and can 
be described with a classical order parameter which corresponds to the wave- 
function of the Bose-condensed atoms. In the noninteracting case all the 
atoms are described by the same single-particle wavefunction which obeys 
the Schrodinger equation 

in^ = -—V^^ + V{r)<l/. (6) 
dt 2m 



In the presence of disorder, this wavefunction can be Anderson localized. Bose- 
Einstein condensates represent an appealing system where it is possible to di- 
rectly study the effect of localization. By using the techniques described in the 
previous section, one is able to create disordered potentials in an extremely 
controlled way, knowing precisely the kind and amount of disorder. Further- 
more, the wavefunction (more precisely, the squared modulus of it) can be 
directly observed by imaging the condensate with a CCD camera. As a result, 
the typical exponential tails of Anderson localized states could be observed (at 
least in principle, if the imaging resolution and sensitivity are good enough), 
allowing the detection of localization. 

An extremely interesting, and still open, problem regards the effect of interac- 
tions on localization. Originally, Anderson formulated his theory for noninter- 
acting quantum particles. If one considers real interacting particles, however, 
the scenario could be significantly different. In the case of electrons, repulsive 
long-range interactions are present due to the Coulomb electric force. In the 
case of ultracold neutral atoms long-range dipolar interactions can be nearly 
always neglected and the dominant interaction mechanism is represented by 
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noninteracting weakly interacting (U<J) strongly interacting (U>J) interactions 

Fig. 5. Pictorial representation of the different interaction regimes for a bosonic gas 
in the presence of disorder. 



s-wave elastic collisions. These short-range interactions can be either attrac- 
tive or repulsive, although only repulsive interactions a llow for the existenc e 



of stable BECs with arbitrarily large number of atoms (jPalfovo et al.l . Il999l ) 



When the BEC density is sufficiently small, as it happens in many experi- 
mental situations, the effect of the coUisional forces can be described within 
a mean-field approach by adding a nonlinear term in Eg. (6), which b ecomes 



the well-known Gross-Pitaevskii equation (GPE) (jPalfovo et al.l . Il999l ): 



K 



ih— = V^^ + \/(r)^ + c/|^r^ . 



dt 



2m 



(7) 



The interaction strength, described hj g = Anh'^a/m, is parametrized as a 
function of one single scalar parameter a, which takes the name of scatter- 
ing length. A similar kind of cubic nonlinearity is present also in the Maxwell 
equations describing the propagation of light in a nonlinear optical medium 
where the index of refraction depends on the light intensity (Kerr effect). This 
term is respo nsible for ma n y effe cts of nonlinear dyn amics, such as soli tonic 



propagation flBureer et a. 



199911. four- wave mixing (jPeng et al.l . Il999l ) and 



instabilities (IWu and Niul . l200ll : iFallani et al.l . |2004J ). The presence of inter- 
actions can heavily affect the physics of localization, which is intrinsically a 
single-particle effect, holding for linear waves. From a naive point of view, 
negative nonlinearities (a < 0, arising from attractive interactions between 
particles, or self-focusing behavior of the wave) could play in favor of local- 
ization. On the contrary, positive nonlinearities (a > 0, induced by repulsive 
interactions, or self-defocusing behavior) are expected to play against localiza- 
tion, making the problem much more interesting to study, both theoretically 
and experimentally. 

The interplay between disorder and interactions in the physics of localization 
has been the object of a very intense theoretical investigation. It was soon 
realized that repulsive interactions can compete with disorder and eventually 
destroy the localization. In strongly interacting systems, however, different 
regimes can be achieved and new quantum phases can be reached in which 
interactions and disorder co-operate in localizing the system in glassy states 
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Fig. 6. Ground state of a BEG in a harmonic + disordered potential. Absorption 
images of the atomic density distribution following a time-of-flight after the re- 
lease from the confining potential. The numbers in the bottom indicate the average 
speckle height Vs in units of the BEG chemical potential ^ ~ 1 kHz in the harmonic 
trap. Adapted from iLve et al.l ^200^ ). 



(see Sec. 4 for a discussion of these phases). An oversimplified picture of the 
different regimes for an interacting bosonic gas in a disordered potential is 
sketched in Fig. 5. Real phase diagrams, of course, are much more complicated 
than this pictorial representation, the details depending e.g. on the kind of 
disorder and on the dimensionality of the system. Actually, there are still many 
open questions to be solved, which the experimental study of interacting Bose- 
Einstein condensates in disordered potentials could address. As a matter of 
fact, these systems offer the advantage of a broad tunability of the Hamiltonian 
parameters, including the kind and amount of the disorder (as seen in the 
previous section) and the interaction strength between the atoms (e.g. by 
using Feshbach resonances fllnouve et al.l . Il998h l 



3.1 A Bose-Einstein condensate in a disordered potential 



3.1.1 Static properties 

A natural starting point to gather information on the behavior of the BEC 
in the disordered potential is the shape of the atomic density distribution af- 
ter release from the confining potential. This time-of-fiight det ection technique 



has b een used since the first experimental realization of BEC (j Anderson et al. 



20051 ) as a precious tool to study its ground state properties. We start con- 
sidering the case of disor dered pote i itials created with optical speckles, first 
investigated with '^^Rb in lLye et al.l (120051 ). 



Basically, one can observe three different regimes, depending on the ratio be- 
tween the speckle height Vs and the BEC chemical potential fi. For very small 
optical potentials < 0.1/i one does not observe any significant deviation 
from the ordinary Thomas-Fermi shape of the BEC expanding from the har- 
monic trap (fig. 6A,B). For higher speckle heights O.l/z ^ V5 < /i one observes 
that the density distribution is strongly modified by the appearance of complex 
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Fig. 7. Ground state of a BEC in a harmonic + disordered potential. In-situ images 
of the atomic density distribution. The sp eckle potential V/? is increasing from top 
to bottom: 0, O.S^u, 0.5/i, 1.3^. Taken from lChen et al.l ^2Q04 ). 



structures in the form of elongated stripes (fig. 6C). Finally, further increasing 
the speckle height to Vs > yU, the expanded density profile ceases to be char- 
acterized by stripes and one can detect only a broad unstructured gaussian 
distribution (Fig. 6D,E). 

The appearance of density modulations in the regime of weak disorder has 
recently attracted a large interest. The structures observed after expansion 
could arise either from real in-trap density modulations or from phase fluc- 
tuations converted in to density modula t ions after time - of-flig ht. More recent 



experimental works (jChen et al.l . l2008l : IClement et al.l . |2008| ) systematically 



investigated such structures, measuring the fringes visibility as a function of 
the speckle height and evidencing that the fringes pattern is stable for the 
same realization of dis order. This observatio n, together with the comparison 
with GPE simulations (I Clement et al.l . l2008l ). suggests that for the actual ex- 
perimental parameters the most plausible scenario is the one in which small 
in-trap density fluctuation s are amplified during the time-of-flight. The prob- 
lem is studied in detail in IClement et al.l (l2008l ) with a thorough analysis of 
the mechanisms involved during the BEC expansion. 

Further increasing the intensity of the disordered potential to V5 > /i the 
condensate is spht up in many condensates localized in the randomly-spaced 
minima of the speckle potential. The absence of interference structures in the 
observed density distribution is due to the fact that the spacing between differ- 
ent condensates is not uniform and gives rise to an interference pattern that, 
averaged over the optical resolution of th e system, is almost flat. This frag- 
mentation scenario has been confirmed in I Chen et al.l (l2008l ) by using direct 
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in-situ imaging of the trapped ^Li atoms (see Fig. 7). 



3.1.2 Collective excitations 



Low-energy collective excitations of the BEC have been studied in the pres- 
ence of a weak speckle disorder producing a corrugation of the harmonic trap 
potential. Since the first production of BECs, frequency measurements of col- 
lective modes have provided precious information for the identification of su- 
pe rfiuidity and, mo re in general, for the characterization of quantum fluids. 
In iLye et al.l (120051 ) the frequency and the damping rate of the dipole and 
quadrupole modes of a BEC in a speckle potential have been measured as a 
function of the strength of disorder. With the term of "dipole mode" one usu- 
ally refers to rigid center-of-mass oscillations of the BEC in the parabolic trap, 
while the "quadrupole mode" for a highly anisotropic trap (or "axial breath- 
ing mode") indicates a shape oscillation in which the center-of-mass does not 
move, b ut the size of the B EC along its l o ng ax is is periodically changing 
in time (jPalfovo et al.l . Il999l ). In iLye et al.l (120051 ) a damping of both modes 
has been observed and measured as a function of the strength of the speckle 
potential. 



Damp ing of the dipole mode has been recently observed also in IChen et al. 
(120081 ) . in which the study of the different dynamic regimes have been ac- 
companied by in- situ detection of the atomic density distribution. Indeed, 
Chen et al.l (120081 ) demonstrated that the halting of the center-of-mass mo- 
tion for strong disorder has to be connected with the creation of a fragmented 
BEC. In this regime different condensates are trapped in the different speckle 
potential wells and no global phase coherence is present, due to the extremely 
long tunneling times between different fragments. 



Lye et al.l (120051 ) measured also a frequency shift of the quadrupole mode. The 
frequency of the quadrupole mode is particularly important since its value 
depends not only on the strength of the parabolic potential but also on the 
Bose-Einstein equation of state, hence it depends on the nature of t he system. 



wheth er it is superfluid or not. The frequency changes measured in lLye et al 



(120051), however, just reflected the change in the effective potential curvature 
induced by the corrugation produce d by the di s order ed potential. The problem 
has been addressed theoretically in iModugnol (120061 ) by numerical solution of 
the Gross-Pitaevskii equation combined with a sum-rules approach, confirming 
the frequency shift as an effect due to the change in the effective trap frequency. 



The effect of the speckle potential produced in these experiments is mostly 
classical, and it does not really produce a change in the nature of the quantum 
fluid. We shall discuss more about this point in the next section, evidencing 
how the correlation length of the potential plays a crucial role for the observa- 
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tion of truly disordered-induced localization effects. In particular, concerning 
the measurements of collective excitations, the presence of disorder with short 
correlation length can modify the superfluid e quation of s t ate le ading to non- 
trivial frequency shifts, as recently studied in iFalco et al.l (120071 ) . 



3.2 The quest for Anderson localization 



In this section we will discuss the state-of-the-art of the experiments aiming to 
observe Anderson-like localization for Bose-Einstein condensates propagating 
in disordered optical potentials. 



3.2.1 Localization in a speckle potential 



Out of the condensed matter systems for which it has been originally proposed, 
Anderson localization has been widely sea rched, and eventually demonstrated, 



in cla ssical wave propagation experiments (IWiersma et al.l . 119971 : ISchwartz et al. 



20071 ). In this kind of experiments, an electromagnetic wave undergoes multiple 
scattering from the randomly-distributed scatterers of the disordered medium. 
Strong (Anderson) localization sets in when the multiple scattered waves in- 
terfere destructively in the propagation direction and localized sta t es be come 
populated. According to the loffe-Regel criterion (lloffe and Regell . Il960l ) this 
happens when the mean free path of the wave becomes as small as its wave- 
length. 

Experiments performed in 2005 at LENS (Florence) and in the group of A. 
Aspect at Institut d'Optique (Orsay) aimed to realize such scattering con- 
figuration with Bose-Einstein condensates propagating in disordered optical 



poten tials produced with speckle patterns (jClement et al.l . l2005l : iFort et al. 



20051 ). The idea behind these two works was quite similar: an initially trapped 
Bose-Einstein condensate of ^^Rb was lef t free t o expand in a one-dimensional 
disordered waveguide. In IClement et al.l (120051 ) this waveguide w as produced 
by a highly elongated magnetic trap, while in iFort et al.l (120051 ) by a single 
beam optical trap. The propagation of the condensed matter wave in the 
waveguide was studied as a function of the height of the disordered speckle 
potential. 



In Fig. 8, taken from iFort et al. 
densate 



(120051 ) ■ the density distribution of the con- 
imaged in situ after a fixed expansion time in the optical waveg- 
uide, is shown for different speckle potential heights (ranging from = to 
Vs = 0.7 fi, with fi the BEC chemical potential) together with the picture of 
the actual speckle field used. Without speckles the condensate freely expands, 
while in the presence of the speckles both the expansion and the center-of- 
mass motion (induced by a small acceleration along the waveguide) start to 
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Fig. 8. Expansion of a BEC in a disordered optical guide. Top) Intensity profile of 
the speckle field used in the experiment. Bottom) Density profiles of the condensate 
after expansion in the disordered optical guide for different speckle heights Vs, here 
expressed in u nits of the BEC chemical potential /i = 2.5 kHz in the initial trap. 
Adapted from lFort et al.l 



be suppressed for Vs > 0.3yU. A closer look shows that actually two different 
components can be distinguished: while a low density cloud expands without 
stopping, a few localized density pea ks become obse r vable when increasing the 
speckle height. In Fig. 9, taken from I Clement et al.l (120051 ). the rms size of the 
BEC expanding in a disordered magnetic waveguide is plotted as a function 
of time for different heights of the disordered potential: one can clearly see the 
transition from a diffusive regime in the absence of disorder to a "localization" 
regime when disorder is present. 



Further investigations have demonstrated that this suppressed expansion is 
not Anderson localization, but a classical localization that can be explained 
with simple energetic arguments. The expanding condensate is not a monochro- 
matic flux of atoms all moving with the same velocity: since the momentum 
distribution of the sample has a finite width (mainly caused by the atom-atom 
repulsive interactions which initially drive the expansion), a low velocity com- 
ponent of the cloud is always present and get trapped in the speckles since it 
has not sufficient energy t o escape the deep est potential wells (in the case of 
red-detuned speckles, as in lFort et al.l (120051 )) or to tunnel thro ugh the highest 



potential barriers (in the case of blue-detuned speckles, as in IClement et al. 
(l2005h ). 
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Fig. 9. Expansion of a BEC in a disordered magnetic waveguide. Time evolution of 
the axial BEC rms size L and center-of-mass position (inset) for various amplitudes 
7 of the random speckle p otential (in units of the BEC chemical potential). Taken 



from lClement et all toO^ ). 



Several theoretical works have studied the e xpansion of an intera c ting: Bose- 



Einstein conde nsate in a speckle potential (Clement et al.l. l2006l: iModugno 



20061: iShapird. 120 07: Sanchez-Palencia et al. 



Sanchez-Palencia et al.l . |2008| ) . In particular, in IClement et al.l (120051 . |2006| ) 



20071: lAkkermans et al.L 



2008 



and Modugn ol (l2006h it has been shown that the expanded BEC density profile 
is actually made up of two spatially separated parts. In the center of the cloud 
interaction energy is dominating over kinetic energy and the BEC density 
profile exactly follows the shape of the potential, as expected from the Thomas- 
Fermi approximation for an interacting Bose gas: no Anderson localization 
is expected to appear in this region. In the wings of the cloud the density 
is much smaller and kinetic energy is dominating over interaction energy: 
here the BEC almost behaves as a noninteracting gas and the density profile 
has deviations from the Thomas-Fermi approximation. However, numerical 
studies based on the Gross-Pitaevskii equation evidenced that no Anderson 
localization is present even in this region. 



There are two possible physical reasons impeding the observation of Anderson 
localization. The first is indeed the presence of interactions: from an intu- 
itive point of view, repulsive interactions between the atoms force them to 
spread more in space, contrasting localization. The second reason is the finite 
correlation length of the disorder: even in the absence of disorder, Anderson 
localization could not be observable because the disordered potential is not 
"good" enough to produce the scattering strength which is necessary to have 
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a localization length smaller than the system size. It is well known that, in 
the pure ra ndom case, any infinite simal amount of disorder leads to localiza- 
tion in ID ( Abrahams et all 1979 ). but in finite-sized systems (like a trapped 
Bose-Einstein condensate) the localization length plays an important role. 



The further experiments reported of iFort et al.l (120051 ) have evidenced that, 
apart from the problem of interactions, the typical speckle potentials employed 
so far in the experiments were not "fine-grained" enough to produce quantum 
refiection/transmission, which is at the basis of ID Anderson localization. 
This has been observed by studying the collision of a BEC with a potential 
defect created with a tightly focused laser beam, that mimics the effects of 
one single speckle grain. The absence of quantum refiection from the potential 
well created with this optical (red-detuned) defect indicated that, in order to 
have quantum scattering, one should use much steeper potentials, i.e. speckle 
potentials in which the autocorrelation length a is smal ler (t he typical cor 



relatio n length of the speckles used in iFort et al.l (120051 ) and IClement et al 
(120051 ) was 5 /im). 



This problem has been theoretically addressed in iModugnd (120061 ). where the 
coefficients of quantum refiection by a potential well and of quantum transmis- 
sion from a potential barrier have been calculated as a function of the potential 
steepness and of the velocity of the incident matter waves. Of course, quan- 
tum refiection/transmission is fundamental in the case of ID localization, in 
which scattering just happens along a line: hence, for multiple scattering to 
appear, only a fraction of the incident wave has to be refiected/transmitted. 
In higher dimensions interference due to multiple scattering could happen also 
for classical refiection from potential hills, for which the requirements are less 
stringent. However, in higher dimensions localization itself is more difficult to 
achieve, owing to the larger localization lengths (if localization is present). 

Despite the obstacles discussed above, Anderson localization could be even- 
tually observed in ID diff usion exper i ments similar to the ones reported in 



Clement et al.l (120051 ) and iFort et al.l (120051 ) . provided that the spec k le au 



tocorrelation length is made small enough. In ISanchez-Palencia et al.l (120071 ) 
the density profile of the BEC expanding in a weak speckle potential has 
been analytically worked out, evidencing exponential localization in the dilute 
tails of the wavefunction, where density is very low and interactions can be 
neglected. In the case of interacting ^^Rb BECs this could happen for very 
small disor der correlat ion l engt h , one order of n i agnitu de less than the ones 
achieved in IClement et al.l (120051 ) and iFort et al.l (120051 ) . The crucial parame- 
ter, as we shall see in the following section, is the ratio between the disorder 
correlation length and the heal ing length, which is the l ength scale associated 
to the effect of interactions. In ISanchez-Palencia et al.l (120071 ) a mobility edge 
was also found as a maximum wavevector of the expanding BEC above which 
localization cannot be observed, which is a peculiar characteristic of speckle 
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potentials with finite correlation length. 



3.2.2 Localization in a bichromatic lattice 



Quasiperiodic lattices, introduced in Sec. 2.2, are a particular class of po- 
tentials which exhibit prope rties common to both periodic and disordered 
systems (jPiener et al.l . 1200 ll ). As in the case of periodic lattices, their spec- 
trum show reminiscence of energy bands. On the other hand, owing to the 
lack of any translational invariance, they support the existence of localized 
states, w hich behave very simi larly to the ones supported by truly disordered 
systems ( iGrempel et al.l . Il982l ). Therefore they can be used as a tool to study 
quantum localization, as a valid alternative to speckle patterns, with the ex- 
perimental advantage of an effortless production of short length-scale potential 
fluctuations. 



Localization in incommensurate bichromatic potentials is a well known topic. 
This problem has been studied in detail in the framework of the Harper model 
(lHarperl . ll955l ) and of the ID tight-binding Aubry- Andre model (lAubry and Andrei . 
1980l ). which is described by the Hamiltonian 



H = -JJ2i\n){n + l\ + |n + + A^cos(27r/?n)|n)(n| 



(8) 



where J is the tunneling rate between next-neighboring sites and A is the 
amplitude of the quasiperiodic modulation of the potential energy, being (3 
an irrational number. The Aubry-Andre model can be experimentally real- 
ized when the primary lattice height Si is much larger than the secondary 
lattice height S2 (see Sec. 2.2). The primary lattice discretizes the system and 
produces a renormalization of the effective mass m* = mi?Ri/J7r^, while the 
height of the secondary lattice S2 = A/Epa is the control parameter which 
drives the localization transition. As a matter of fact, differently from what 
happens with a pure disordered potential, in the quasiperiodic local- 
ization transition exists even in ID, with a critical value A 2J of the 
quasidisorder amplitude for producing a localized ground state. 

This behavior is illustrated in Fig. 10a, in which we plot the lowest-energy 
single-particle eigenstate in the incommensurate lattice, as obtained by nu- 
merical integration of the ID Schrodinger equation (6), which holds in the 
continuum. The figure shows, in grayscale, the squared modulus of the ground 
state wavefunction as a function of position (horizontal axis) and disordering 
lattice strength (vertical axis). One clearly sees that for low values of disorder 
the ground state is an extended state (see Fig. 10c), i.e. the wavefunction ex- 
tends across the entire lattice (in the limit A = one recovers the extended 
Bloch state describing the ground state in a periodic lattice). Increasing dis- 
order above a threshold value the wavefunction suddenly localizes around few 
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position [\xm] 

Fig. 10. Transition from extended to localized states in a bichromatic incommensu- 
rate potential, a) The square modulus of the ground state wavefunction is plotted in 
grayscale as a function of position and disordering lattice strength. b,c) Logarithmic 
plots of the ground state wavefunction below and above the localization transition. 



lattice sites, with exponentially decreasing tails, in the same way as in the case 
of Anderson localized states in 5-correlated disordered potentials (Fig. 10b). 

So far, we have shown that the incommensurate lattice supports the existence 
of localized states for noninteracting particles. What hap pens when one intro- 



duces interactions? This problem has been discussed in iSchulte et al.l (12005 
20061 ). where the ground state of the system has been calculated for a three- 



colour optical lattice. Introducing repulsive interactions between the atoms, 
the numerical integration of the ID Gross-Pitaevskii equation shows that the 
ground state wavefunction becomes a superposition of many single-particle 
localized states, which add up to form an overall extended state, as shown 
in Fig. 11 for different interac tion streng t hs. Si milar results for a bichromatic 
lattice have been presented in lLye et al.l (120071 ). 



This behav ior can be interpreted in terms of a screening effect induced by in- 
teractions ( Sanchez-Palencia . 20061 : Schulte et al. . 2006). The nonlinear term 
in the Gross-Pitaevskii equation (7) can be treated as an effective potential 
cancelling the spectral components of the original potential varying on length 
scales larger than the healing length C, = 1/ ^/Siian. This latter quantity is the 
typical length scale that is asso ciated to the variatio n of the BEG wavefunction 
around sharp potential jumps (jPalfovo et al.l . Il999l ). More generally, the heal- 
ing length is the typical length scale below which the condensate wavefunction 
is able to behave quantum-mechanically. As a consequence, in order to observe 
localization, the interference effects producing localization should take place 
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Fig. 11. Ground state of the Gross-Pitaevskii equation for a Bose-Einstein con- 
densate in a three-colour optical lattice for different effective interaction strength 
g = 0.2 (a), ^ = 8 (b), ^ = 256 (c). The latter value corresponds to the interac- 
tion strength for a Bose-condensed sample of = 10^ ^^Rb atoms in an elongated 
harmonic trap with frequencies lo±_ = 27r x 40 Hz and u)\\ = 2tt x 4 Hz. Taken from 
Schulte et aiT(|2005l ^. 



on a distance smaller than the healing length, otherwise the BEC wavefunc- 
tion would behave classically over longer distancies. The healing length can 
be made larger by reducing the amount of interactions in the system, that 
could be achieved either by reducing the scattering length a or by reducing 
the atom density n. If the healing length is smaller than the disorder local- 
ization length < (, see Eq. 5) no localized states can be observed. If the 
healing length is larger than the localization length, but smaller than the sys- 
tem size (C < ^ < L), one could observe a superposition of many localized 
states. Finally, if the healing length is the largest length scale in the system 
> L > () the BEC wavefunction collapses in a single localized state. This 
crossover is illustrated in Fig. 11, where the ground state in the three-color 
lattice is plotted for the same lattice heights but different interaction strength 
(hence same ( but different ^): while in panel a) > C and a few localized 
states can be easily detected with clearly exponentially decreasing tails, in 
panel c) ^ < ( and an overall extended state forms. 



The existence of localized states can be probed with t ransport experiiri ents, 
similar to those presented in the previous section. In iLye et al.l (120071 ) the 
transport of an ^^Rb BEC has been studied in the presence of a bichromatic 
incommensurate potential. Localization of the center-of-mass motion has been 
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Fig. 12. Oscillations of a BEC in a parabolic + incommensurate two-color lattice, a) 
Center-of-mass position for N = 1.5 x 10^ atoms, si = 10 and different disordering 
lattice heights S2 = 0.1 (filled black circles), S2 = 0.17 (filled gray circles), S2 = 0.25 
(empty circles), b) Center-of-mass position for si = 10, S2 = 0.25 and different atom 
num bers N = 1 -5 x 1 0^ (empty circles) and N = 2 x 10^ (filled circles). Adapted 
from lLve etaP (|2007l ). 



observed, the stronger the smaller is the strength of interactions (tuned by 
changing the atomic density), as shown in Fig. 12. This density-dependent 
behavior, with interactions pushing to delocalize the system, is reminiscent of 
Anderson-like localization. However, in the regime of parameters studied in 
this work, no simple Anderson-like localization has to be expected, being the 
eigenstates of the system similar to the state dominated by interactions shown 
in Fig. 11c. The suppression of the center-of-mass motion shown in Fig. 12 
was mainly caused by the strong modulation of the BEC wavefunction on the 
length scale of the beating between the two lattice periods, which resulted in 
very low tunneling times across the lattice and, consequently, in an extremely 
slow dynamics. 



Effects of nonlinear dynamics have also been considered as possible mech- 
anisms to damp the motion. As a matter of fact, suppression of transport 
is expected to appear for r nechanisms alike the in t eraction-indu c ed dy nam- 
ical instability observ ed in ICataliotti et al.l (120031 ) , iFallani et al.l (120041 ) and 



Cristiani et al.l (120041 ) for monochromatic optical lattices. In the system stud- 



ied in these works the interplay between repulsive nonlinearities and band 
structure resulted in fast-growing excitations dephasing the system and halt- 
ing the motion. The same effect can also be observed in bichromatic optica l 
lattices, in which a band structure can still be identified (jPiener et al.l . l200ll ). 
with a multitude of energy gaps opening in the spectrum and getting denser 
and denser with increasing height of the secondary incommensurate lattice. 



From what discussed above we can draw a preliminary conclusion. Bichro- 
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matic potentials do allo w to solve t he pr o blem of the co r relation length of 



the sp e ckles produced in iLye et al.l (120051 ): IClement et al.l (120051 ): iFort et al. 



( I2OO5I ): ISchulte et al.l (120051 ) . which causes the localization length to be too 
large to be observable. In bichromatic potentials the localization length can 
easily be smaller, however the presence of too strong interactions still remains 
and makes it impossible to observe a clear localization of the wavefunction 
in a few localized states. In order to achieve clear signatures of Anderson 
localization one needs to work with extremely weakly interacting samples. 



3.3 Further directions 



From the theoretical point of view, as we have already pointed out, the effect 
of disorde r on the interacting B ose gas is an extremely interesting topic of re- 
search. In iLugan et al.l (l2007al ) the ground state of an interacting Bose gas in 
a disordered potential has been deeply studied. In particular, a Lifshitz glass 
phase has been introduced characterizing the ground state of the system for 
weak interactions. Lifshitz states are a particular class of localized states which 
exhibit "weaker" localization properties than Anderson-localized states, in the 
sense that they show exponential decay only in the very far tails, while close to 
the maximum their shape mostly depends on the local properties of the poten- 
tial. This means that they mostly resemble bound states of isolated potential 
wells or trapped states between barriers, differently from Anderson- localized 
states, whose shape is determined by global properties of the potential, i.e. by 
the combined effect of many impurities / potential wells. In Fig. 13 we show 



the pha se diagram of the interacting disordered BEC derived in iLugan et al. 



(j2007al ) as a function of the BEC chemical potential and of the speckle height. 
Starting from this Lifshitz glass phase and increasing interactions, a phase of 
fragmented interacting BECs has been proposed, which is a precursor of the 
Bose glass phase (see Sec. 4 for further discussion). 



The BEC fragmented state has been previously described in lWang et al.l (120041 ) 
In this work the authors studied the ground state of a BEC in the disordered 
potential produced by the random imperfections of a magnetic microtrap. This 



pape r was motivated by several experirnental observations (IFortagh et al. 



2OO2I : iLeanhardt et all . l2003l : lEsteve et all . l2004l : Ijones et all . 120041 ) , in which 



fragmentation of the BEC at very close distances from t he current-carrying 



wires of the microchip was observed. In the same work (jWang et al.l . |200J) 



the BEC dynamics in the disordered potential was also investigated. In par- 
ticular, the spectral analysis of the sloshing motion after displacement of the 
confining potential allowed to identify different dynamical regimes: superfiuid 
oscillations, self-trapping and an intermediate chaotic regime. 



Recent theoretical works (IBilas and Pavlofa . l2006l : iLugan et al.l . l2007bl ) have 
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Fig. 13. Phase diagram of an interacting disordered BEC as a function of the chem- 
i cal po tential /x and of the disordered potential height Vr. Taken from lLugan et al 



also studied the problem of Anderson localization of excitations in a Bose- 
Einstein condensate. In the weakly interacting case, BE C excitations are de- 



scribed by the Bogoliubov theory (jPalfovo et al.l . Il999l ). By calculating the 



leading-order many-body corrections to the classical BEC wavefunction, one 
finds that, in the absence of external potentials, excitations are described by 

'{hcky + {h^ky2my, with k 



quasiparticles with dispersion relation hj = 
wavevector of the excitation and c sound velocity inside the BEC. This spec- 
trum has two distinct regions with different ^-dependencies. For k -C (with 
^ the healing length) the excitation spectrum is phonon-like and excitations 
have energy uj ^ ck. For k ^ the spectrum is particle-like and the energy 
of the excitations is a; ^ Uo + hk^ /2'm, w here Uq = mc^/h is an energy shift due 



to interactions. 



Bilas and Pavlofll (120061 ) have shown that, in the presence of a 



white-noise random potential, excitations can undergo Anderson localization, 
almost in the same way as the whole BEC wavefunction can und ergo. Further 



detailed studies have been carried out in iLugan et al.l (l2007bl ) in the more 
realistic case of correlated disorder, showing that the localization length (and, 
correspondingly, the possibility to observe localization in finite-sized BECs) 
crucially depends on the correlation length of the disorder. 



The effect of disorder on the coherent BEC dynamics can be observed also on 
the dephasing of Bloch oscillations which is expected to appear when a disor- 
dered or quasi-disordered pote r itial i s superimposed on a tilted optical lattice 
( ISanchez-Palencia and Santod . l2005l : ISchulte et al.l . 120081 ). Bloch oscillations 
are the coherent oscillations of a wavepacket in a periodic potential when a 



25 



constant force is applied. This phenomenon has been observed for the first 
time with ultracold atoms in optical lattices (IRaizen et al.l . 119971 ) because of 
the much longer coherence times than the ones achievable for electrons moving 
in real solid-state lattices, where defects and impurities strongly dephase the 
system in a time much shorter t han the oscillation p eriod. Also interactions 
lead to dephasing, as evidenced in lMorsch et al.l (200l|) when Bloch oscillations 
were observed for the first time in a Bose-Einstein condensate. Later it was 
demonstrated that this interaction-induced dephasing can be controlled and 
eventually cancelled by tuning the interaction strength wi th F eshbach reso- 
nance s, as recently demonstrated in lGustavsson et al. fl2008l) andlFattori et al 



( I2OO8I ) , or by using ultracold fermionic samples (IRoati et al.l . l2004j ) , for which 
interactions are forbidden by the Pauli principle. Noninteracting particles in 
perfectly periodic optical lattices perform undamped Bloch oscillations and, 
thanks to this possibility, they can be used as micr oscopic probes fo r high - 



precision measurements of forces at small distances (ICarusotto et al.l . l2005l ). 



Starting from this ideal situation and adding disorder on top of the periodic 
lattice, one can quantitatively study the deph asing induced by dis order in a 
controlled way, as first theoretically st udied in ISchulte et al.l (120081) and then 



recently investigated experimentally in iDrenkelforth et al.l (|2008l ). 



Generally speaking, disorder leads to the disruption of coherent effects. Quite 
interestingly, however, under certain conditions disorder can induce a sponta- 
neous ordering of the system. This effect, known as random-field induced or- 
der, has been originally studied in the context of classical spin models, which 
in the presen ce of disorder may exhibit a magnetiz ation higher than in the 



ordered case (IWehr et al.l . l2006l : ISen De et al.l . 120071 ). Two-component Bose- 



Einstein condensates in the presence of a random Raman coupling between 
th e two states can be used t o study this class of effects, as recently proposed 
in 



Niederberger et al.l (120081 ). 



3.4 Anderson localization: the state of the art 



Starting from the first experiments w i th speckle potential s ( iLye et al.l . 12005 



Fort et all . l2005l : IClement et al.l . l2005l : ISchulte et al.l . 120051 ) , the quest for An 



derson localization in Bose-Einstein condensates has been a strongly active 
direction of research. On one side experimental groups have focused o r i pro- 



ducing di s order ed potentials on thinner length scales (IClement et al.l . 12006 



DeMarcol . l2007l ). in order to increase the amplitudes of quantum scattering 



and decrease the attainable localization lengths. On the other side, the chal- 
lenge is to reduce atom-atom interactions in order to make localization ob- 
servable. Once this is obtained, it will be even more interesting to study the 
effect of adding a controlled amount of interactions. For this purpose it can be 
strongly helpful to take advantage of Feshbach resonances to tune the scatter- 
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ing length a, which is the key parameter defining the strength of interactions 
g = Anh'^a/m (see Sec. 3). In this perspective, the choice of the element un- 
der investigation is c rucial. The first experiments performed with BECs in 



disordered potentia ls (ILye et al.l . l2005l : iFort et al.l . l2005l : I Clement et al.l . 12005 



Schulte et al.l . l2005l ) have focused on Rb, which is a quite convenient element 



for the implementation of cooling schemes, but has the disadvantage of hav- 
ing a quite large scattering length a ~ lOOao (with oq the Boh r radius) and 



no fa vorable Feshbach resonances at convenient magnetic fields f lMarte et al. 



2OO2I ) . A much easier tuning of atom-atom interactions could be provided by 



different elements, such as ^Li (which has b een already studied in combina- 
tion w ith laser speckles in I Chen et al. dioosh ) or K, studied by iRoati et al.l 
(l2007h . 



While this review was being completed (march 2008) two experiments have 
succeeded in observing Anderson localization of coherent matter waves, in the 
groups of A . Aspe ct in France and here at LENS. In the French experiment 
(IBilly et al.l . |2008| ) a Bose-Einstein condensate is left free to expand in a dis- 
ordered waveguide produced by combining a weakly focused laser beam with 
a ID speckle potentia l. Di fferently from the conceptually similar experiments 



Clement et al.l (120051 ) and iFort et al.l (120051 ) , the analysis of the in-situ den- 



sity profiles shows clear indication of exponentially decreasing tails, which is 
a signature of Anderson localization. This have been made possible by a com- 
bination of several factors: the small atomic density in the tails (necessary to 
reduce the counteracting effect of interactions), the small speckle autocorre- 
lation length (necessary to have many quantum scattering events during the 
diffusion) and the high detection sensitivity (allowing the observation of expo- 
nential decay of the density). The obser ved localization is then quanti t ativel y 
compared with the theory developed in ISanchez-Palencia et al.l (120071 . l2008l ). 



In the experiment at LENS (IRoati et al.l . |2008| ). Anderson localization has 



been observed for a noninteracting K BEC in an incommensurate bichro- 
matic lattice, similar to that used in lLye et al.l (120071 ). Here the strategy to ex- 
clude the effect of interactions is different: instead of working with dilute sam- 
ples, interactions are cancelled by tuning a static magnetic field in proximity 
of a Feshbach resonance to set the scattering length to zero. The noninteract- 
ing condensate in the quasiperiodic potential thus realizes the noninteracting 
tight-binding Aubry- Andre model of Eq. (8), which exhibits a transition from 
extended to localized states for increasing disorder. The crossover between 
extended to localized states is studied in detail by looking at the expansion of 
the BEC and by studying spatial and momentum distribution of the states, 
all of which result in agreement with the Aubry- Andre predictions. 
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4 Strongly interacting regime 



In the previous section we have discussed the physics of disordered weakly- 
interacting bosonic systems. The theoretical description of this regime is pro- 
vided by the semiclassical Gross-Pitaevskii Eq. (7), which describes the prop- 
agation of nonlinear matter waves. When interactions are strong, however, 
this mean-field description is not capable to fully explain the behavior of the 
system. A more appropriate description is provided by a full quantum the- 
ory, taking into account quantum correlations between particles. Also in this 
strongly interacting regime disorder may induce localized quantum phases: 
these have a different nature from Anderson localization since correlations be- 
tween particles are important, whereas Anderson localization is essentially a 
single-particle effect. 

Experimentally, a convenient way to enter the strongly interacting regime is 
provided by the use of optical lattices (which we have already introduced 
in Sec. 2.2). In a deep optical lattice the system becomes effectively stronger- 
interacting because of the combined effect of the tighter squeezing of the atom 
wavefunction in the potential wells (with a consequent increase of the local 
density) and of the increase in the effective mass due to the finite tunnel- 
ing times across the potential barriers (which makes the kinetic energy less 
important with respect to the interaction energy). 

For a system defined on a lattice, starting from the full many-body Hamilto- 
nian, one can derive a simplified zero-temperature model, in the approximation 
that all the particles occupy the fundamental vibrational state of the lattice 
sites. In this limit the quantum state of an interacting gas of identical bosons in 
a lattice pote ntial is well describ e d by the second qu antization Bose-Hubbard 



Hamiltonian (IFisher at al.l . Il989l : iJaksch et al.l . Il998l ) 



h = -jY. + ^ E % - 1) + E (9) 

{hi') 3 j 



where bj (6^) is the annihilation (creation) operator of one particle in the j- 
th site, fij = b^jbj is the number operator, and (j, j') indicates the sum on 
nearest neighbors. Each of the three terms on the right-hand- side of Eq. (9) 
accounts for a different contribution to the total energy of the system: J is 
the hopping energy, proportional to the probability of quantum tunneling of 
a boson between neighboring sites, U is the on-site interaction energy, arising 
from atom-atom on-site short-range interactions (repulsive for ^''Rb, for which 
U > 0) and giving a nonzero contribution only if more than one particle occu- 
pies the same site, and ej G [—A/2, A/2] is a site-dependent energy accounting 
for inhomogeneous external potentials superimposed on the lattice. 
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Fig. 14. Qualitative phase diagram for a disordered system of lattice interacting 
bosons. Three phases can be identified: a superfluid (SF), a Mott insulator (MI) 
and a Bose glass (BG). 



The quantum phase of the system depends on the interplay between these 
three energy scales: hopping energy J, interaction energy U and disorder A. 
We start considering the ideal case of a translationally invariant system, in 
which A = 0. Assuming integer filling of the sites, when J > U the system 
is in a superfluid (SF) state, in which the bosons are delocalized across the 
lattice and the tunneling ensures off-diagonal long-range coherence. Instead, 
when U > J, the system is in a localized Mott insulator (MI) state, where 
long-range phase coherence is lost and number Fock states are created at the 
lattice sites. The actual phase diagram of the system depends on the chemical 
potential (related to the atomic density ) and shows the ex istence of MI lobes 
with integer number of atoms per site (iFisher at al.l . Il989l ). In the left graph 
of Fig. 14 we show a qualitative sketch of the phase diagram for a 3D system. 



The transition from a SF to a MI for ult racold bosons in an optical lattice 
has been proposed i n iJaksch et al. fll998h and reported for the first time in 
reiner et al.l (|2002| ). where the ratio J/U was varied across the transition 
point by controlling the height of the lattice. The realization of a MI state does 
require a 3D optical lattice, since, in order to enter the strongly interacting 
regime, the atomic wavefunction should be squeezed in tightly confining traps, 
with a site occupation on the order of unity. However, by using deep optical 
lattices effectively slicing the atomic sample into decoupled 2D or ID systems, 
it is possible to s tudy the S F-MI transition injower dimensionality, as made 



m 



Stoferle et al.l (120041 ) and lSpielman et all (120071 1. 



In the presence of a disordered external potential the additional energy scale 
A enters the description of the system and is responsible for the existence of 
a new quantum phase. In the presence of weak disorder the MI lobes in the 
phase diagram should progressively shrink and a new Bose glass (BG) phase 
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should appear (central graph of Fig. 14) , eventually washin g away the MI 
region for A > U (right graph of Fig. 14) (IFisher at al.l . 119891 ). In a simplified 
view, a Bose glass is half-way from a Mott insulator to a superfiuid: it is an 
insulating state, with no long-range phase coherence, as the Mott insulator 
is; nevertheless, it is compressible and has no energy gap in the excitation 
spectrum, as a superfiuid has. 



The Bose glass phase has been first identified in lGiamarchi and Schulzl ( 1l988l ). 
where strongly interacting ID bosonic systems were studied. In the '90s it was 
widely studied in the context of the superfiuid-insulator transition observed 
i n many condensed-m atter systems, such as ^He ad sorbed on porous media 



(ICrowell et al.l . Il995l ). thin supercond ucting films (IGoldrnan and Markovic 



19981 ). arrays of Josephson ju nctions (Van de r Zant et al. . 199^ ) and high- 



temperature superconductors (jjiang at al.l . 1 1994 : Budhani et al.l . Il994j ). The 



possible realization of a Bose glass in a system of ul t racol d bosons in a disor- 



dered lattice has been first proposed in lDamski et al.l (120031 ) : lRoth and Burnett 



( 120031 ). More recently, the phase diagram of this system has been derived in 



other theoretical paper s , consid ering also finite temp erature effects ( iKrutitskv et al 



2006 : Buonsante et al. . 2007a ). detection schemes ( Bar-Gill et al. . 20061 ) and 



the p ossible realizati o n of a Bose gla s s with incommensurate bichromatic lat- 



tices (jPugatch et al.l . l2006l : iRoscildd . 120071 ). Evidences for Bose glass phases 



have been also theoretically obtained for different classes of B ose-Hubbard 



model s, where disorder is introduced either i n the hopping energy (IBuonsante et al. 



2007bl ) or in the on-site interaction energy (iGimperlein et al.l . boosh 



The Bose glass is just the simplest disordered quantum phase that can be 
realized in the strongly interacting regime. When atoms of different species, 
or different internal (spin) states of the same species, are considered, more 
complicated models can be experimentally realized and new disordered quan- 
tum phases can emerge. Atomic Bose/Fermi mixtures, in particular, represent 
a ver satile system in whi c h many differen t disor dered models can be real- 



ized ( ISanpera et al.l . l2004l : lAhufinger et al.l . |2005| ). In the strong interacting 



limit this system can be described in terms of composite fermionic particles 
corresponding to on e fermion + one bosonic particle/hole in the same site. 



Sanpera et al.l (120041 ) have shown that the interaction between these compos- 



ite fermions can be tuned by changing the external potential: thus, a disordered 
potential can be used to induce an effective random interaction between the 
particles. This possibility allows the investigation of a variety of disordered- 
related r nodels, from fermionic Ising spin glasses to models of quantum per- 
colation (lAhufinger et al.l . 120051 ) . 



30 




Fig. 15. Excitations in the deep insulating phases, a) In a Mott insulator the tunnel- 
ing of one boson from a site to a neighboring one has an energy cost A.E = U. b) In 
the disordered case the excitation energy is AE = [/ it Aj , that becomes a function 
of the position. In the Bose glass state, in which |Aj| > U, an infinite system could 
be excited at arbitrarily small energies and the energy gap would disappear. 

The quest for Bose glass 

Experiments with disordered bosons in the strongly interacting regime started 
at LENS in 2006. The system under investigation was a collection of ID atomic 
systems in a bichromatic optical lattice. A main optical lattice was used to 
induce the transition from a weakly interacting superfluid to a strongly cor- 
related Mott insulator. A secondary optical lattice was then used to add con- 
trolled quasi- disorder to the perfect crystalline structure of the MI phase. 
With reference to Eq.( 9), the non-commensurate periodic potential superim- 
posed on the main lattice introduces inhomogeneities of the energy landscape 
ej G [—A/2, A/2] on the same length scale as the lattice spacing. 



4.1.1 Excitation spectrum and coherence properties 

As introduced in the previous section, the excitation spectrum is an important 
observable that can be measured in order to characterize the quantum state 
of the system. By exploiting the possibilit y of ti me-modulating the lattice 



potential, as first realized in IStoferle et al.l (12004J ). it is possible to directly 



measure the excitation spectrum and study how it is modified by the presence 
of disorder. Naively speaking, in the MI phase one realizes a crystal of atoms 
pinned at the lattice sites and sitting on the fundamental vibrational level, 
as schematically shown in the top of Fig. 15. In a MI an energy gap in the 
excitation spectrum exists, since the elementary excitation - the hopping of a 
particle from a site to a neighboring one, or, in other words, the creation of 
a particle-hole pair - has an energy cost U, corresponding to the interaction 
energy of a pair of mutually repelling atoms sitting on the same site (see Fig. 
15). 



31 



In Fig. 16a we show the excitation sp e ctrum of a Mott insulator measured in 
the LENS experiments (IFallani et aLl . 120071 ) . The plot shows a well resolved 
resonance at energy U, which is distinctive of the MI state, and a second reso- 
nance at energy 2U. While the physical origin of the excitation peak at U is the 
tunneling of particles between sites with the same occupancy, the second peak 
at 2U can be ascribed to several processes: it can arise from tunneling at the 
boundary between MI regions with different site occupancy (that are present 
due to the inhomogeneity of the confined sample), from higher-order processes 
and from nonlinear effects due to the strong modulation. A theoretical anal- 
ysis of the response of the bosonic system to t his lattice modula t ion ha s been 
recently reported in lKoUath et al.l (120061 ) and IClark and JakschI ( l2006l ). 



When increasing disorder the experiment showed a broadening of the reso- 
nance peaks, which eventually become undistinguishable when A ^ U. As a 
matter of fact, the presence of disorder introduces random energy differences 
Aj G [—A, A] between neighboring sites (see bottom of Fig. 15). As a con- 
sequence, the tunneling of a boson thro ugh a potential barrier costs U ± Aj, 
that becomes a function of the position ( Guarrera et al. . 20071 ). The excitation 
energy is not the same for all the bosons, differently from the pure MI case, 
and the resonances become inhomogeneously broadened, as can be observed 
i n the experimental spectra at weak disorder (A < U) shown in Fig. 16b, c 



(IFallani et al.l. l2007f). This br oadening is in agreement with a semi-classical 



cal works (IHild et al.l . 12006 



model (iGuarrera et al.l. 120071) and has been recently predicted in theoreti 



20081 ). where the authors study the 



IZakrzewskil . 

dynamical response of a ID bosonic gas in a superlattice potential when a 
periodic amplitude modulation of the lattice is applied. 



Eventually, when A > U, one expects that an infinite system can be excited at 
arbitrarily small energies and that the energy gap would shrink to zero. When 
this happens, nearby sites become degenerate and regions of local superfluidity 
with short-range coherence appear in the system. This novel many-body state 
in which there is no gap but the system remains globally insulating is a Bose 
glass. 



From the experimental point of view, additional information on the nature of 
the many-body ground state can be acquired by analyzing the density distri- 
bution of the atoms released from the lattice after a time-of-fiight. Long-range 
coherence in the sample results in a density distributio n with interferenc e 



peaks at a distance proportional to the lattice wavevector (jPedri et al.l . l200ll ) 



The visibility of these peaks provides a measurement of phase coherence. When 
increasing the height of th e main lattice, a prog ressive loss of long-range coher- 
ence has been reported in IFallani et al.l (120071 ) indicating the transition from 
a superfiuid to an insulating state, also in the presence of disorder. The com- 
bination of the excitation spectra measurements and the time-of -fiight images 



indicates that, with increasing disorder, the system realized in IFallani et al. 
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Fig. 16. Excitation spectra of the atomic system in a Mott insulator state for in- 
creasing height of the disordering latti ce. The resonances a re lost and the excitation 
spectrum becomes flat. Adapted from lFallani et al.l (j2007l ). 



( 120071 ) goes from a MI to a state with vanishing long range coherence and a 
flat density of excitations. The concurrence of these two properties cannot be 
found in either a SF or an ordered MI, and is consistent with the formation 
of a Bose glass, which is indeed expected to appear for A > U. 



Much work has still to be done for the exhaustive characterization of such novel 
disordered state. New detection schemes should be implemented, in order to 
have access to additional observables. This necessity is not only restricted 
to the study of disordered systems, being a more general issue shared by 
the experimental investigation of different strongly interacting lattice systems, 
including e.g. systems with magnetic orderii ig or mix t ures of different species . 
From the theoretical side, very recent works (jRoscildd . 120071 : IRoux et al.l . l2008l ) 
have extensively studied the problem of ID interacting bosons in quasiperiodic 
lattices, working out the phase diagrams (which include the presence of Bose 
glass and incommensurate "band insulating" / "charge density wave" regions) 
and studying how the different phases affect experimentally detectable signals. 



4- 1-2 Noise correlations 



In the recent work iGuarrera et al.l (120081 ) noise interferometry has been used 
to study interacting ^^Rb bosons in the bichro r natic lattice. This detection 
technique, originally proposed in lAltman et al.l (120041 ). is based on the anal- 
ysis of the spatial density- density correlations of the atomic shot noise after 
time-o f-flight. These correlations are ba sed on the Hanbury Brown & Twiss 
effect (jHanbury Brown and Twissl . ll956l ): if two identical particles are released 
from two lattice sites, the joint probability of detecting them in two separate 
positions (e.g. imaging them on two separate pixels of a CCD camera) depends 
on the distance between the detection points. These correlations, arising from 
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quantum interference between differe nt detection paths, were first observed 
for bosons in a Mott ins ulator state ( Foiling et al.l . |2005| ) and then also for 
band- insulating fermions (IRom et al.l . l2006l ). The sign of the correlations de- 
pends on the quantum statistics: while bosons show positive correlations (due 
to their tendency to bunch, i.e. to arrive together at the detectors), fermions 
exhibit negative correlations (due to the antibunching, consequence of the 
Pauli exclusion principle). In the case of a bosonic Mott insulator, one ob- 
serves positive density-de nsity correlation pea ks at a distance proportional to 
the lattice wavevector ki flFoUing et al.l . |2005| ). as shown in the first image of 
the bottom row of Fig. 17 for the recent experiment at LENS. 



In lGuarrera et al.l (120081 ) noise correlations have been measured, starting from 
a Mott insulator state, for increasing heights S2 of the secondary lattice. The 
absorption images after time-of-flight do not present significative differences, 
as shown in the top row of Fig. 17, and demonstrate the absence of first order 
(phase) coherence of the atomic system in the insulating state, even in the 
presence of the secondary lattice. However, second order (density) correlations 
turn out to be significantly different with varying S2, as illustrated in the 
noise correlation functions plotted in the bottom row. More precisely, with 
increasing S2, one observes the appearance of additional correlation peaks at a 
distance proportional to the wavevector k2 of the secondary lattice and to the 
beating between the two lattices ki — k2- These peaks have to be associated 
with the redistribution of atoms in the lattice sites as the disordering lattice 
is strengthened: the MI regions characterized by uniform filling are destroyed 
and atoms rearrange in the lattice giving rise to a state with non-uniform 
site occupation, which follows the periodicity of the secondary lattice. The 
redistribution of atoms is then quantitatively detected by measuring the height 
of the additional correlation peaks. 



Noise correlations thus prove to be a tool to extract important information on 
the lattice site occupation, which is connected to the second-order correlation 
function of the many-body state. The appearance of si milar correlation p eaks 
was predicted in t heoretica l works for hard-core bosons (IRey et al.l . 120061 ) and 



soft-core bosons (iRoscildd . 120071 ) in bichromatic lattices. Future works will 
study the possibility to use noise interferometry to get additional insight on 
the nature of the disordered insulating states produced in the experiment, in 
particular in connection with the realization of a Bose glass phase. 



4-2 Experiments with atomic mixtures 



As we have discussed in Sec. 2.3, disorder can be produced by letting the atoms 
interact with randomly-di stributed scatterers of a di fferent atomic species. The 
configuration proposed in iGavish and CastinI (120051 ) has not yet been realized 
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Fig. 17. Top) Time-of- flight absorption images of atoms in a Mott insulator state for 
increasing height of the secondary lattice S2- Bottom) Density-density correlation 
functions corresponding to the pictures above. The additional correlation peaks for 
large ,52 arise from the destruction of t he Mott domai n s and the redistribution of 
the atoms in the lattice. Adapted from Guarrera et al. ( 20081 ) . 



experimentally. However, in 2006 the first experiments with binary mixtures 
in optical lattices hav e been realized alm ost at the same ti me in two differ- 



ent g roups, in Zurich (iGiinter et al.l . l2006l ) and in Hamburg (lOspelkaus et al. 



20061 ). In these experiments ^^Rb bosons and fermions were mixed together 
in a 3D optical lattice. Since the two atoms have very similar resonance wave- 
lengths, the depth of the optical lattice is almost identical for the two species, 
however, being potassium lighter than rubidium, its mobility is favored. As 
a result, for a range of lattice heights, rubidium can be localized in a Mott 
insulator state, while potassium atoms are still able to move across the lattice. 

In these experiments the superfluid to Mott insulator transition of '^^Rb was 
investigated as a function of the concentration of impurities, typically in 
the range to 20%. In particular, the visibility of the interference pattern after 
time-of-fiight was investigated. The observation reported by the two groups 
was a downshift of the lattice height value at which coherence starts to be 
lost, when potassium atoms are introduced in the system. Different interpre- 
tations for this effect have been given, including finite-temperature effects, 
disorder-like induced localization, or effects connected with the strong attrac- 
tive interaction between the two species. 



Recently, a closely related system has been investigated at LENS (ICatani et al. 
20071 ) by using a binary bosonic mixture of ^^Rb and ^^K in a 3D optical lattice 



Similarly to the experiments described above, loss of coherence in the rubid- 
ium sample induced by the presence of potassium has been observed. However, 
this latter experiment differs from the former ones in two points: the mixture 
is bosonic/bosonic (instead of bosonic/fermionic) and the interspecies inter- 
action is repulsive (instead of attractive). The observation of similar effects in 
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systems with different quantum statistics and different interaction sign rufes 
out some of tlie interpretation given so far, even if a clear explanation of the 
observation has not yet been found. Future advances of these experiments with 
mixtures in optical lattices will be pushed by the use of Feshbach resonances 
for fine tuning of the interspecies interaction. 



5 Conclusions 



The investigation of Bose-Einstein condensates in disordered potentials is a 
fastly growing field of research. For decades condensed-matter physicists have 
theoretically studied the interplay between disorder and interactions in deter- 
mining the transition from metals to insulators. Now disordered systems can 
be realized in cold atoms laboratories and, differently from traditional solid 
state systems, they allow a fine tuning of both disorder and interactions, as 
well as the advantage of new detection capabilities, thus extending the range 
of the possible experimental investigations. 

The first experiments were realized only a few years ago, with the successful 
creation of disordered and quasi-disordered potentials and the first studies of 
the behavior of ensembles of ultracold bosons in different regimes of inter- 
actions. Although this field of research is quite young, it already relies on an 
extensive literature, mostly comprising theoretical works. The physics of disor- 
dered atomic systems is indeed extremely rich, both in the weakly interacting 
regime (where Anderson localization and its disruption by interactions can be 
studied) and in the strongly interacting regime (where particles are strongly 
correlated and new quantum localized phases can emerge, as the Bose glass) . 
The few experiments performed until now have just opened a new direction, 
showing the great potentialities of ultracold atoms for investigating the physics 
of disorder, but still leaving open questions concerning the observed localiza- 
tion effects. Regarding this point, much more results are likely to come in the 
near future. 

We have tried to give an overview of this newborn field of research, discussing 
the topics of interest and the experimental efforts made up to now. We would 
like to conclude by noting that any review is by necessity incomplete and 
cannot be exhaustive of all the work made in the field. For this reason we 
apologize with the authors of works which we have unintentionally forgotten 
to mention. 
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